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Abstract 

In this letter in the framework of a simple see-saw scenario with two (or three) 
quasi degenerate Majorana neutrinos we proposal that one of these neutrinos 
can be very weakly coupled, yet there is a mechanism of the generation of 
the abundance of such "dark" neutrino in the early universe. The mechanism 
of production is due to oscillations between "dark" Majorana neutrino and 
one of the Majorana's which has relatively large Yukawa couplings ("bright" 
Majorana neutrino). The transition of "bright" Majorana into a "dark" one is 
nonadiabatic. We point out on the similarity with the Landau-Zener transition 
regime. In our model one can explain the observed dark matter density, present 
matter-antimatter asymmetry and active neutrino data all at the same time. 



lr rhis letter is based on the original work pQ. 



1 Introduction 



The observed mater-antimatter asymmetry and the presence of yet unknown form of 
matter in the universe known as Dark Matter (DM) both require the extension of the 
Standard Model which does not contain neither a DM candidate nor a mechanism to 
produce large enough baryon asymmetry in order to match the observed value. On the 
other hand, the observed neutrino mass scales and mixings can not be incorporated within 
the Standard Model (SM) as well. The sew-saw mechanism [2j provide probably the most 
natural and simple explanation why neutrinos are massive, yet so light compared to all 
other Standard Model particles. At the same time, neutrino mixing data represent a 
positive test for leptogenesis |3|, an attractive explanation of the the observed matter- 
antimatter asymmetry of the Universe as a direct consequence of the see-saw mechanism. 
Can DM be accommodated within the see-saw model as well? It is known, that it is 
enough to have just two Majorana's to account for both leptogenesis and active neutrino 
data (see, e.g. [1]). Assuming that the number of Majorana neutrinos is the same as 
the number of lepton flavors in Standard Model one can ask the question if one of the 
Majorana neutrino can be made almost decoupled and, thus, stable at present but being 
produced efficiently enough in the early universe. However, if a particle is almost stable 
at present there is typically quite hard to produce it during earlier epoch. One clear 
possibility is the DM production from the inflaton decay along with the rest of particles. 
There is no argument why not to assume, for example, that inflaton couplings to the DM 
particle and to the SM matter are not of the same order. In that case DM particle can be 
produced as efficiently as ordinary matter. However, such mechanism, due to the present 
lack of the detailed picture of inflation and most importantly the exit from it, is strongly 
model dependent. Besides, nothing would require one of the Majorana neutrinos to be a 
DM particle. Any other particle decoupled from the SM fields would be as plausible. It is 
appealing, therefore, to look for another more constrained mechanism of the production 
of the "dark" Majorana neutrino. The perturbative production out of the thermal bath 
can be a solution. However, typically this would require the reheating temperature after 
inflation to be quite high. 

In this letter we propose that "dark" Majorana neutrino is produced via oscillation 
mechanism which is similar to the Landau-Zener type transition. The production starts 
and terminates near some critical temperature T c due to the oscillations between "dark" 
Majorana neutrino and one of the strongly coupled (or "bright") Majorana neutrino. The 
critical temperature can be as low as the Majorana neutrino mass scale which, in turn, 
can be as low as O(TeV). Contrary to the perturbative production mechanism in our 
model high reheating temperature is not required. In order for oscillation mechanism to 
be successful both neutrinos have to be quasi degenerate in massH We will discuss two 
scenarios. The one, called "minimal", has only the see-saw type couplings being added to 

2 In the case of quasi degenerate neutrinos the Majorana mass scale can be made quite low as the 
consequence of resonant enhancement of the kinematic part of the expression for the CP-asymmetry, 
which compensates for the smallness of the Yukawa couplings [5]. 
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the SM Lagrangian. In this case the production rate turns out to be not efficient enough. 
In the " next-to-minimal" scenario we will introduce a simple dimension five operator which 
couples Higgs field to Major ana neutrinos and is suppressed by the high energy physics 
scale A. Such operator could result from GUT extensions of the SM, supergravity etc. 
This operator is nearly irrelevant for any physics at present. However, at relatively high 
temperatures (T > Tew) it will turn out to be an effective source for the "dark" Majorana 
production H 

2 The Model. 

The see-saw extension of the Standard Model is described by the Lagrangian 

C = C SM + mm - F ik LiN k $ c - ^N?Ni + h.c. (1) 

For simplicity we will assume that only one of the "bright" Majorana neutrino serves as 
a source for the production of the "dark" Majorana neutrino. We will employ the density 
matrix formalism |1Q[ lllj . The evolution of the number densities of the Majorana and 
of the active neutrinos is described by the density matrix which obeys quantum kinetic 
equation: 

ipu = [Uk, Pk] - + r£, P } + irl (2) 

where T^ p are destruction and production rates correspondingly (see, e.g. [12]). In gen- 
eral case, this is 2 x (3 + Nm) matrix equation, where Nm is the number of the Majorana 
neutrinos. In practice, one can always neglect the off-diagonal elements pnl of this matrix 
because the latter undergo rapid oscillations. Because the helicity flip and lepton number 
violation processes are suppressed the opposite helicity parts of this matrix evolve inde- 
pendently. This simplifies the Eq. ([2]) further. Finally, in our model the oscillations of 
"bright" neutrino into a "dark" one is effective only during a brief period, so that light 
lepton part of the density matrix does not have any effect on the final abundance of the 
"dark" Majorana neutrino. Therefore, it is enough to consider only the part of this matrix 
which describes two quasi degenerate Majorana neutrinos: one "dark" and one "bright". 
This is described by 2 x 2 matrix the p = (p)ab is (A, B are flavor indices). This equation, 
in general, has to be solved for each mode k and the resulting spectral matrix p^ then has 
to be integrated over the phase volume to obtain the number density matrix p: 

P = J (2^3 P* (3) 

3 A simultaneous picture of baryogenesis and DM has been already proposed [BJ|7]. Here the matter- 
antimatter asymmetry is generated through RH neutrino oscillations [3] and the lightest O(KeV) mass 
RH neutrino is a warm DM. The other two Majorana neutrinos in this scenario are much heavier but still 
have to be lighter then the electroweak scale as well. Previously, some other models were already proposed 
where the light sterile neutrino is produced by the mixing with the ordinary light ones and plays the role 
of Dark Matter [jj . In this letter we explore the case when Majorana neutrinos mass scale is well above 
the electroweak scale. 
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The diagonal elements of the matrix p give the physical number densities of the two 
neutrinos after being multiplied by the massless fermion number density (or, equivalently, 
by the photon number density n~). 

In our model one of the neutrinos is decoupled and the other is in thermal equilibrium]! 
In that case the Eq. ([2]) simplifies to 

ip\ = [-ffk,Pk]- (4) 

Furthermore, instead of solving the Eq. ([3]) for each mode we will adopt a more simplified 
approach choosing k ~ T in order to facilitate further discussion. A more quantitative 
analysis will be given elsewhere [Tj. In this letter we aim on a qualitative discussion so 
that the latter simplification is justified. Because the abundance of the "dark" Majorana 
neutrino is expected to be much smaller then thermal abundance of the "bright" Majorana 
neutrino at the moment of production, we can set pbb = P°bb- ^ n addition we will assume 
that there is no CP violation in the A — B sector. With all these assumptions the Eq. (J3|) 
for a given momentum k ~ T leads to the following set of equations: 

PAA = 2H A bP-, 

P- = P+{Hbb ~ H A a) ~ H A bP C bb' 

P+ = P-{Haa~ H B b), 

where we have defined p + = (pab + Pba)/2, P- = —i (pba — PAb)/2- The Hamiltonian 
can be computed from the real part of the Majorana neutrino two-point Green function 
on Fig. (1). Below we will consider two cases with two different effective temperature 
dependent Hamiltonians result from self-energy corrections shown on Fig. (2-3). 

3 The Minimal Scenario. 

The Lagrangian of the model ([!]) is written in the Majorana neutrino mass eigenstate 
basis. The Yukawa matrix in the case of two generations is given bj0 

(F) { =( Fn Fl2 )^( ° Fl2 ) (5) 
13 V F22 ) v F21 F22 J ' 

where by assumption Fn <C Fi.2,21 *C -^22, therefore we neglect Fn keeping only Fi2,2i 
elements. The Hamiltonian can be computed following [13] and is given b)H 




(H) ij = jM? + k2I + - Km K , 8 AM|, (6) 



4 We assume that the critical temperature T c satisfies the condition M < T c < T cq , where M is the 
Majorana neutrino mass scale and T cq is the temperature at which the second neutrino come to the thermal 
equilibrium. This condition put some constraints on the model, that will be discussed later. 

5 The indices z,j will be used for the mass eigenstate basis, while letter indices will represent flavor 
basis. The index A is always reserved for the "dark" Majorana neutrino. 

6 Assuming that the two RH neutrinos are strongly degenerate, such that A = \M% — M\\/M\ 1. 
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Figure 1. Majorana neutrino two-point Croons function. 
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Figure 2. Correction to the Majorana neutrino self-energy: minimal scenario. 
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Figure 3. Correction to the Majorana neutrino self-energy: next-to- minimal scenario. 
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where I is the identity matrix and the matrix K is defined as 



K = F ] F. (7) 

The first term, proportional to the identity matrix, cancels out in the evolution equation 
so that we will not keep it further. The computation of the "dark" Majorana neutrino 
production can be made in the mass eigenstate basis by computing pn and noting that 
Pli ~ PAA- Instead, we will make a rotation into the flavor basis and compute paa directly. 
The diagonalization of the matrix Fij can be made with the help of two unitary matrices 
U (rotation of the Majorana neutrino basis) and V (rotation of the light lepton basis). 
Introducing two different angles Q\j and 9y and diagonalizing matrix F^ one obtains 

{F)ab*>( J. > (8) 



F 22 

where it was used that both angles Q\jy are small and given by 

Ouy - ^ (9) 

The Hamiltonian is rotated according to 

H' = U^HU, (10) 

and one obtains the Hamiltonian familiar from the active neutrino oscillations and MSW 
effect: 

^'W-(^ _M^| |rA . 2 J. (») 

where E = yjM 2 + k 2 , M = Mi ~ M 2 . Assuming that E ^k^T, AM 2 2 X = -2M 2 A < 0, 
where A = (Mi — M2VM2 <C 1 one obtains solving the Eqs. ([5]) 

PAA *W ISA) Kll -2M Kn > (12) 



where m p i = M p i y^90/8ir 3 g^ ~ 10 18 GeV. The last inequality follows from the requirement 
that the critical temperature T c = M(8A/i^22) 1//2 a t which H' BB becomes equal to zero 
satisfies T c > M. 

It is instructive to compare (]12p to the value of the Landau-Zener parameter Q that 
controls the regime in which the transition occurs. Its value is given by (see, e.g. [14]) 

M 2 A sin 2 26 u 

Q ~ Eco^OuhR ' (13) 
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where 



hR 



2T 



2T 2 
mpi 



(14) 



For MSW effect [15J it is essential that Q > 1 while in the Landau-Zener regime Q < 1. 
Taking into account account that ^ « 1 and estimating for k ~ T the value of (fl3l) one 
obtains at T = T r 



Q 



rn.pi 
M 



K 



22 
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K 



XX- 



(15) 



One can see that pn ~ Q, as expected since both represent the probability of the transition 
in the nonadiabatic regime. The upper bound on the DM energy density can be estimated 
as follows 



10 25 ^ii. 



0.1 



(16) 



Given that the lifetime of the flavor state Na has to exceed the age of the universe in 
order to be a DM candidate one can estimate that K u < 1(T 42 (1 TeV/M x ) and the 
corresponding density (|16|) is by far too small. 



4 The next-to-minimal scenario. 

As it was shown in the previous section the minimal model fails to produce the correct 
amount of the "dark" Majorana neutrino abundance. In this section we propose an ex- 
tension to the previous setup by introducing (explicitly) no new particles but allowing a 
dimension five effective operator 

C cS = { -^\^N c A N B , (17) 

where is usual Higgs field, A is some high energy scale and the elements of the ma- 
trix (X)ab ~ 0(1)- The correction of the type (|17|) is generic if one assumes a certain 
additional physics at energies much higher then the electroweak scale. For example the 
operator (|17p generically will result from supergravity with A being M p i or from tadpole 
type graphs in GUT theories via heavy Higgs exchange. There are other possibilities to 
generate (]17p . In what follows we will not address any specific model which provides the 
relation of the scale A to any other scales like M p i , Mqut , -^susy etc. Instead we will 
treat A as phenomenological parameter and impose the constraints on its value from the 
requirement that the operator (]17p is responsible for the correct DM abundance at present. 
Additionally, we assume that 

< 18 > 
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The latter conditions simply implies that the the operator (|17p in the Majorana neutrino 
flavor basis gives a contribution to the off-diagonal entrees of the Hamiltonian which 
dominate over the term 

AM 2 sin 29 , x 

— 51— (19) 

at all temperature above the electroweak phase transition. In that case we may disregard 
(fl9l) in the computation below. The validity (fTSl) can be easily inforced by choosing Ojj 
small enough (or even setting it to zero). 

The resummed Majorana propagator will contain the nondiagonal entrees which are 
given by 

- S AB « 1 ^S(q 2 - ml)n h (q), (20) 

where 

MQ) = , (21) 

e t — 1 

with u being a four velocity of the thermal bath. The resulting expression for the real 
finite part of T*ab in Boltzmann approximation is 

T 2 

- RepAB - Z AB (T = 0)] « ^-ZiKiMup-f, (22) 

where z$ = (m$/T) ~ 0(1), and i^i(^) is a modified Bessel function. Alternatively, we 
can neglect thermal mass of the Higgs. In that case one obtains 

T 2 

- Re[£ AB - Yiab{T = 0)] « Ym^ 7 "' (23) 

In what follows the exact coefficient in front of T 2 /A, which is of order of O(10~^) is not 
particularly important and we will use (|23l) for simplicity. Assuming as well thatHI 

\k 22 » ^. (24) 



results in the following effective Hamiltonian: 







ji2 



" * i T 2 M 2 A ! 2 | T ^ 22 • (25) 



12A £ k 



Solving Eq. @ one obtains that 



m pl M (oa s 

P AA A2i^' ( 26 ) 

7 This assumption is not crucial but simplifies analytics and turns out to be satisfied for a wide range 
of parameters after restrictions coming from the "dark" neutrino abundance at present and its decay rate 
are imposed. 
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where z c is given by 

r ^22 



8A 

The above estimate transfers into the DM abundance which is at present 

^DMh 2 1 /10 24 GeV\V M \ 



0.1 " z c V A J Vl0 9 GeVy 



(27) 



(28) 



The operator (|17|) after electroweak symmetry breaking leads to the mixing term 
(v^ w /A)NaNb between "dark" and "bright" neutrinos. This results in the decay of 
the "dark" neutrino via oscillations into the "bright" one and its consequent decay. The 
rate of this process is 



T{N A ^N B ^ ...) = S _ (29) 



The decay rate of the "bright" neutrino is dominated by T(Nb — > $ + and one estimate 
that 

TiNA ^ ...) ~ (™) ! (^)\v. (30) 

There is phase space of parameters {M, A, z c } which leads to the correct DM abundance 
and the lifetime of the "dark" neutrino larger then the age of the Universe. For example, 
one can choose z c = 10 -2 , M = 10 7 GeV and A = 10 24 GeV. For this choice of parameters 
the condition (|24|) is satisfied as well. The "extreme" case would be to take z c = 10 -3 
(because, at smaller values it is hard to have the "bright" neutrino being thermally pro- 
duced), M = 1 TeV and A ~ 10 22 5 GeV. In this case the "dark" neutrino has the mass 
near electroweak scale and a lifetime close to the present age of the Universe. Of course, 
a more stringent bounds may come from astrophysical observations and the lower bound 
on the neutrino mass scale quoted here may go up substantially. 

There is an another constraint which comes from the four-body decay Na — ► 2$ + 
Nb — > 3$ + v. The rate of such decay can be estimated as 

7 M 3 24 / M \ 4 /i0 24 GeV\ 2 

F 4 ~ UT' W K* ~ 10- 24 (^^j ^^—j eV. (31) 

Requiring that this rate is smaller then the value of the Hubble parameter at present leads 
to the constraint 

This implies that M can not be higher then about 10 8 — 10 9 GeV depending on the value 
of Zn. 
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5 Conclusions 



In this letter we proposed the see-saw model in which one of three Majorana neutrinos 
("dark" neutrino) is weakly coupled and is a DM candidate. The mechanism of the 
generation of a suitable abundance was suggested. At least one of the strongly coupled 
Majorana neutrino ("bright" neutrino) have to be quasi degenerate with the "dark" one. 
The allowed quasi degenerate Majoranas mass scale M can, in principle, be as low as 
O(TeV) and up to as high as usually suggested in GUT scenarios, i.e. M ~ 10 s — 10 9 GeV. 
The mass of the third neutrino is adjusted correspondingly: if M is below 10 8 GeV the 
third neutrino has to be also quasi degenerate with the other two in order to have lepton 
number produced efficiently Otherwise it remains a relatively free parameter as in the 
standard hierarchical leptogenesis scenario. This letter is based on the results which will 
be reported in much more details in [T]. 
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